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Abstract 

A commonly employed measure of the signal amplification properties of an input/output 
system is its induced C? norm, sometimes also known as H^o gain. In general, however, it 
is extremely difficult to compute the numerical value for this norm, or even to check that it 
is finite, unless the system being studied is linear. This paper describes a class of systems 
for which it is possible to reduce this computation to that of finding the norm of an associ- 
ated linear system. In contrast to linearization approaches, a precise value, not an estimate, 
is obtained for the full nonlinear model. The class of systems that we study arose from the 
modeling of certain biological intracellular signaling cascades, but the results should be of 
wider applicability. 



1 Introduction 

The analysis of signaling networks constitutes one of the central questions in systems biology. 
There is a pressing need for powerful mathematical tools to help understand and conceptualize 
their information processing and dynamic properties. One natural question is that of quantify- 
ing the amount of "signal amplification" in such a network, meaning in some sense the ratio 
between the size of a response or output and that of the input that gave rise to it. See for 
instance [5 1 for a recent paper in this line of work. 

In control theory, a routine way to quantify amplification is by means of the induced C? 
norm or "-ffoo gain" of a system. A major difficulty when trying to apply these techniques to 
signaling networks is that such systems are usually highly nonlinear. Thus, typically, mathe- 
matical results are only given for small inputs or "weakly activated" systems, see for instance 
lI5i|3J. For large signals, that is, when analyzing the full nonlinear system, even deciding if the 
norm is finite or not is usually a very hard question. 

In this paper, motivated by the particular systems studied in ['5l,'3l, we introduce a class of 
nonlinear systems, which includes all these motivational examples as well as many others, and 
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we show finiteness and how to obtain precise values for norms, by reducing the problem of 
norm estimation to the same problem for an associated linear system. This associated system 
is sometimes a linearization of the original system around an equilibrium point, though it need 
not be. In any case, the techniques are not at all related to linearization techniques, but instead 
borrow from comparison theorems, ISS-like estimates, and the theory of positive systems. 

2 Definitions and Statements of Results 

We deal with systems of the following special form: 

x{t) = A{x{t)) x{t) + B{x{t)) u{t) , x{0) = (1) 

(or just ",x = A{x)x + B{x)u"), where dot indicates time derivative, and states x{t) as well as 
input values u{t) are vectors with nonnegative components: x{t) e M>q and u{t) e M>q for all 
i > 0, for some positive integers n and m. We view A and B as matrix valued functions 

where M>q = (M>o)*^, for any positive integer k, is the set of vectors ^ e M'^ in Euclidean 
/c-space with all coordinates > 0, i = 1, . . . , A;. Associated to these systems we also have an 
output or measurement 

y{t) = h{x{t)) = C{x{t))x{t) 
taking values y{t) e W, for some integer p, where C : R>q — > R^^". 

Assumptions 

We make several assumptions concerning the matrix functions A, B, and C, as follows. 
Stability: 

The matrix A{0) is Hurwitz, that is, all eigenvalues of A{0) have negative real parts. 

Maximization at ^ = 0: 

For each ^ G R^q, A{^) < A{0), B{^) < 5(0), and C(^) < C(0), meaning that 
A{Oij < ^(0)ji for 6^ch i.j e {1, . . . B{^)ij < B{0)ij for each i e {1, . . . ,n} and 
j G {1, . . . ,m}. andC(,^)jj < C{0)ij for each i e {1, . . . ,p} and j e {1, . . . ,m}. 

Positivity of system: 

For each ^ G M>o and each i G {1, . . . , n} such that = 0, it holds that: ^(Oij ^ foi" 
all j ^ i and B{^)i~> for all j. Also, for every ^ G Qj(^) > for all i, j. 

Local Lipschitz assumption: 

The matrix functions A{^), B{^), and C(^) are locally Lipschitz in ^. 
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Remarks about the form of the system 

The special form assumed for the system is in itself not very restrictive, since every (affine 
in controls) system x = F(x) + B(x)u may be written in this fashion, provided only that 
F be a continuously differentiable vector field and F(0) = 0, for instance by taking A(^) = 
F'{X^) dX, where F' indicates the Jacobian of F. This reduction to a "state dependent linear 
form" X = A{x)x + B(x)u is often useful in control theory, where it appears for instance in 
the context of "state-dependent Riccati equation" approaches to optimal control. Of course, 
the difficulty is in satisfying the above assumptions for A and B. 

A special case in which these hypotheses are satisfied is that of models of cell signaling 
cascades as in |l5l|3l. These are systems whose equations can be written as follows (with n 
arbitrary and m = 1): 

Xi = aiu(ci — Xi) — PiXi 

Xi (y.iXi—\(yCi Xi^ Pi-^i 5 ^ 2, . . . , ?1 

and output y = Xn, and the a^'s. Pi's, and q's are all positive constants. We represent this 
system in the above form using: 

° \ 
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Note that A{i) < A{0) and B{^) < 5(0), for all ^ G RIq, because -aS < for all i. The 
matrix A(0) is lower triangular with negative diagonals, and hence is Hurwitz. Positivity holds 
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as well: if ? = 1 and ^ is such that = 0, then = for all j 7^ 1 and B{^)ii = 

(X\C\ > 0; if instead i > 1 and ^ is such that = 0, then = for all j ^ {i — 1, i}, 

= '^jCj > 0, and B(^)ii = 0. Finally, the functions A(-) and B(-) are linear, and 
hence Lipschitz. The matrix C(^) = (0, 0, . . . , 0, 1)^ is constant and nonnegative. Thus all 
properties hold for this example. 

A linear one-dimensional system Xn+i = x„ — ixn+i may be cascaded at the end, as in [l3l|, 
and the output is in that case redefined as y = x„+i; this may be again modeled in the same 
way, and the assumptions still hold. 

Induced gains 

Assume given a system dU). We consider the operator T that assigns the solution function x to 
each input u. To be more precise, we consider inputs u E /^^([O, 00), M>q), and define x = Tu 
as the unique solution of the initial value problem ([T]). In principle, this solution is only defined 
on some maximal interval [0, T), where T > depends on u\ however, we will show below 
that T = +00, and that x is again square integrable (and nonnegative), so we may view x as 
an element of >C^([0, 00), M"q) and T as an (nonlinear) operator 

T : /:2([o,oo),R^o)^^'([0,oo),My . 

We will write |-| for Euclidean norm, and use ||-|| to denote norm: = \u'f dt. For 
the operator T, we consider the usual induced operator norm: 

||T|| := sup . . . 

We will show that ||T|| < 00 for the systems that we are considering. In order to see this, we 
first consider the linear system 

i = A(0)z + S(0)m, z(0) = (2) 

with output V = i{z) = C{0)z, and its associated operator 

L : C\[0, 00), M^o) C\[0, 00), M%) : u ^ z . 

Since A(0) is a Hurwitz matrix, z{t) is defined for all t > 0, and L indeed maps into C^. 
Furthermore, its induced norm the "i^oo gain" of the system with output y = z, is finite; 
see for instance ||H. (The i^oo gain is defined for arbitrary-valued inputs u E jC^([0, 00), R™); 
we will remark below, cf . Section |5J that the same norm is obtained when only nonnegative 
inputs are used in the maximization.) Moreover, the C°° (or induced gain is also 

finite. Therefore, using ||-||^ to denote supremum norm \\z\\^ = supi>o \z{t)\, we can pick a 
common constant c > such that 

< c||m|| and < c for all m g £^([0, 00), M^o) (3) 

where c upper bounds both and \\L\\^ (we use \\L\\^ for operators to denote induced 
£2 ^ £oo ^orm). 
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Our object of study are the compositions with the output maps, i.e. the input/output opera- 
tors: 

To : /:2([o,oo),R^o)^^'([0,oo),R|o) 
: u ^ y = C{x)x = C{Tu)Tu 

and 

: u^v = C{Q)z = C{<d)Lu 
and their corresponding induced norms. Our main result is as follows: 

Theorem 1 The norm of To is finite, and \\To\\ = \\Lo\\. 

3 Preliminary Results 

We start our proof by remarking that the solutions of dU) remain in R>q. To see this, we need 
to verify the following property (this is a standard invariance fact; see for instance [2] for a 
discussion in a related context): 

for each i = 1, . . . , n, each ^ G M>q such that = 0, and each /i G ]R>0' 

Since = 0, we need to prove that ^j^i ^{Oij^j + ^iOijl^j nonnegative, but this is 
implied by the positivity assumption. 

Similarly, solutions of O remain in M|o, as also (A(0)^ + 5(0)/i). > if = 0. 

The next observation is a key one: 

Lemma 3.1 Every solution of ([T]), with u G C^, is defined for all t > 0. Moreover, for any 
two solutions x of ([l]) and © with the same input u, it holds that < Xi{t) < Zi{t) for each 
coordinate i = 1, . . . , n and each t > 0. 

Proof. We use the following comparison principle for differential equations. Suppose that 
f{t,^) and g{t,^) are such that /i(t, ^ Qii^, for alH = 1, . . . , n and all ^ G M>o, and 
that we consider the solutions of i = /(t, x) and z = g(t, z) with the same initial condition 
(or, more generally, initial conditions a;(0) < z(0)). Then, provided that g is quasi-monotone 
(and suitable regularity conditions hold, as here), we may conclude that x{t) < z{t) (compo- 
nentwise) for alH > for which both solutions are defined. See for instance 0161. Quasi- 
monotonicity means that dgi/d^j > for all i j. 

Let us now take any fixed control and let f{t, ^ = ^(0^ + B{C)u{t), g{t, ^ = ^(0)^ + 
B{0)u{t). We have that f{t,^) < g{t,^) coordinatewise, because A{^) < A{0) and B{^) < 
B{0) by assumption. To see that g is quasi-monotone, one needs to verify that A(0)ij > 
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for all i ^ i- but this follows from the positivity assumption on [A, B). Thus the comparison 
principle tells us that x{t) < z{t) for all t > for which the solution x is defined (the solution 
z is defined for all t, since ^ is linear and ^4(0) is a Hurwitz matrix). We already observed 
that X is bounded below by zero; thus, the maximal solution x is bounded on any finite interval, 
and hence it is indeed defined for all t, and the Lemma follows. I 



Corollary 3.2 For each u E C}, the solution Tu of ([T]) is in C}, and the operator T has finite 
norm. Moreover, 

||Tm|| < \Iju\ < c \\u\\ 

and 

\\Tu\\^<\\Lu\\^<c\\u\\ 

where c is any constant as in ©, so in particular ||T|| < ||L|| < c and ||T||^ < \\L\\^ < c. 
Similarly, the i/o operator To also has finite norm, \\Tou\\ < \\Lou\\ and ||ToM||gQ < ||-^^ow||oo for 
all u G C\ and ||To|| < \\Lo\\, \\To\\^ < \\Lo\\^. 

Proof. Pick any u, and let x = Tu and z = Lu. By the Lemma, < Xi{t) < Zi{t) for all t, so 

poo ^ POO ^ 

llxll^ = / 'S^Xi{syds< / Zj^sY ds = \\z\f . 

Jo Jo 



So ||Tn|| < \\Lu\\ < c \\u\\, and since u was arbitrary it follows that ||T|| < Similarly, 

Iklloo = sup|a;(t)| < sup\z{t)\ = \\z\\^ 

t>0 t>0 

leads to II Tm 1 1 < ||Lm|| and ||T|| < ||L|| . 

II lloo — II lloo II lloo — II lloo 

The positivity and the maximization properties for C imply that, for each coordinate i of the 
outputsy(t) = C{x{t))x{t) andv{t) = (:7(0)2(t), wehave < yi{t) = 'E'^=iCij{x{t))xj{t) < 
Yl]=i C'ii(0);Zj(t) = Vi{t), so the inequalities for To and Lo follow by an analogous reasoning. 



Note that the inequality 1 1 To || < ||Lo|| gives the finiteness statement as well as one-half of 
the equality in the main theorem. 

For any matrix Q, we denote by |(5| its induced operator norm as an operator in Euclidean 
space, that is, the smallest constant d such that \Q^\ < d\^\ for all ^. 

Lemma 3.3 There is a nondecreasing and continuous function M : M>o M>o such that: 

\A{o-Am<Mm\^\ 

m)-Bio)\<Mm\^\ 
ic(o-c(o)i<M(iei)iei 

for all C G R^o- 
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Proof. This is a simple consequence of the local Lipschitz property. On each ball B{R) = 
I 1^1 < R}, we pick the smallest common Lipschitz constant Mq{R) for A{-), B{-), and 
C(-). The function Mq is nondecreasing, and hence can be majorized by a continuous and 
nondecreasing function M. Since ^ G we have that \A{^) - A{0)\ < M(|^|) |^|, and 

similarly for B and C. I 

Corollary 3.4 For each function x e C^f] C°°: 

\\A{x{.)) - Am < M{\\x\\J 11x11 

\\Bixi-))-Bm<Mi\\x\\J \\x\\ 
||C(a;(-)) - 5(0)11 <M(||x||J \\x\\ 



where M is as in Lemma 133 
Proof. We have: 



\\A{xi-)) - AiO)f = / \A{xis)) - A{0)f ds 







oo 



< / M{\x{s)\y\x{s)\' ds 
Jo 

POO 

< / M{\\x\\J'\x{s)fds 
Jo 

\x{s)fds 

= M{\\x\\r\\xf 



and similarly for B and C. 



4 Proof of the Main Result 

Pick any input u E and consider once again the respective solutions x = Tu and z = Lu. 
By Corollary 13.21 we know that both < c\\u\\ and < c||m||. Therefore, using 

Corollary 13. 41 we also have that: 



\\A{x{-))-Am<cM{c\\u\\)\\u\\ 

\\B{x{-)) - fi(0)|| < cM{c\\u\\) \\u\\ 
\\C{x{.))-Cm<cM{c\\u\\)\\u\\ 
where M is as in Lemma l331 Let ip : M>o M>o be the function Lp{t) :- 

(A(0) - A{x{t))) x{t) + (5(0) - B{x{t))) u{t) . 
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By the Cauchy-Schwartz inequality, 



||(A(x(.))-A(0))x(-)|| < P(x(.))-A(0)||||x|| 

< c^M{c\\u\\)\\uf 



and 



ii(fi(x(.))-fi(o))oii < mx{.))-Bm\H\ 

< cM{c\\u\\)\\uf 

from which we conclude that 

y\\ < lihw) \\u\\ 

with 7(r) = (c^ + c)M{cr)r, and 7 is a function of class /C, i.e. continuous, strictly increasing, 
and with 7(0) = 0. 

Consider the difference w{t) = z{t) — x{t). Note that w{0) = 0. Evaluating w = [A{0)z + 
B{0)u] — [A{x)x + B{x)u] and rearranging terms, 

w{t) = A{0)w{t) + ip{t) . 

Using once again that A(0) is a Hurwitz matrix, we know that, for some constant d > which 
depends only on A(0) and not on the particular input u being used, || < d \\ip\\. Therefore, 
\\u!\\ < 7(||m||) after redefining 7(r) := d'j(r). 

In terms of the outputs y = ToU = C{x)x and v = LoU = C(0)2;, 

\\v~y\\ = \\C{0)z-C{x{-))x\\ 

< \\Cmz-x)\\ + \\{C{0)-Cix{-))x\\ 

< |C(0)| ||^-x|| + \\CiO)-C{xi-))\\ \\x\\ 

< \CiO)\^i\\u\\)\\t4+c'Mic\\u\\)\\uf 

and we can again write the last term as 7(||tt||) if we redefine 7(r) := |C(0)|7(r) + 
c^M{cr)r. 

The triangle inequality gives us that — \\Tu\\ < \\Lu — Tu\\ and ||i^ow|| — \\Tou\\ < 
\\LoU — Tou\\, and CoroUarv 13.21 gives \\Tu\\ < \\Lu\\ and \\Tou\\ < \\Lou\\, so we may sum- 
marize as follows: 

Proposition 4.1 There is a function 7 G /C such that 

< - ||Tn|| < 7(||u||) 

and 

< 11-^^0^11 — ||Tom|| < 7(||m||) 
for any input u e C^. □ 
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To conclude the proof of Theorem [H we must show that \\To\\ > \\Lo\\. Let g = \\Lo\\, and 
pick a minimizing sequence Un, n = 1, 2, ... of nonzero inputs in that is, 

hm — — — = g . 

n^oo \\Un\\ 

Pick a sequence of real numbers > such that Vn := SnUn — > (for example, e„ = 
[n Since Lo is a linear operator, ||Lot>„|| = ||LoM„||, and since ||f„|| = Sn \\un\\, 

also II Lot'n II / II f nil = II LoMn 11/ II II- Applying the second inequality in Proposition 1431 



< i^-MM < 7(||,;„||)^0 



which gives that ^^T°"n" g, and therefore ||To|| > c, as desired. 



5 Positive vs. arbitrary inputs 

We have shown that the norm of the nonlinear system ([T]) can be exactly computed by finding 
the norm of the associated linear system Q. The computation of induced norms for linear 
systems is a classical area of study, and amounts to the maximization, over the imaginary axis, 
of the largest singular value of the transfer matrix of the system (the Laplace transform of the 
impulse response), the norm; see for instance |4|. There is, however, a potential gap in the 
application of this theory to our problem, namely, the usual definition of Hoo norm corresponds 
to maximization over arbitrary inputs u E >C^([0, oo), M™), not necessarily inputs with values 
in M>Q as considered in this paper. We close this gap now, by showing that the same result 
is obtained, for systems Q, whether one optimizes over arbitrary or over nonnegative inputs. 
We give two proofs, one elementary and the other one less trivial but leading to a stronger 
conclusion. 

The positivity assumptions imply that the operator Lo is a nonnegative convolution opera- 
tor: 

{Lou){t) = [ W{t~ s)u{s)ds, (4) 
Jo 

W{t) e {R>oy''"'yt>o. (5) 

Here W{t) = C(0)e*^'^°)i?(0), and its nonnegativity follows from the fact that e*^ has all 
entries nonnegative, provided that Fij > for all i ^ j. (This last fact is well-known: it is 
clear for small t from the expansion e*^ = I + tF + o{t), and for large t by then writing e*^ 
as a product of matrices e^*/*^-*^ with the positive integer k large enough.) We next show that 
any operator as in (HHSj) has the same norm whether viewed as an operator on >C^([0, 00),]^*") 
or on /^^([O, oo),M>q). Since the norm as an operator on nonnegative inputs is, obviously, 
upper bounded by the norm on arbitrary inputs, it will be enough to show that, for each w E 
£^([0, oo),]R™), there is another input w G £^([0, 00), M>o) with ||w|| = ||'iZ;|| and \\Low\\ < 
\\Low\\. 
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Given such a w, we start by writing w = u — v, where u and v are picked in £^ ( [0, cxd) , M>o) 
and orthogonal. (Such a decomposition is always possible. We define coordinatewise, for 
each i = 1, . . . ,m, Ui := max{wj,0} and Vi := max{— w^.O}; clearly, w = u — The 
supports of Ui and Vi are disjoint, so (wj, Vi) = Ui{t)vi{t) dt = for each i, and also then 
{u, v) = J^iLii^ij ^i) = 0-) We now let w := u + v. Since m and v (or — f ) are orthogonal, 
ll""^!!^ = ll""!!^ + ll""^!!^ = ll""!!^ + ll"^!!^ = ll""^!!^' ll""^!! = ll""^!!- Because Lo is nonnegative, 
both X — LgU and y — LgV are nonnegative. To finish the proof, we only need to see that 
\\x — y\\ < ||x + y||: 



A different proof, which in fact also implies that the supremum in the definition of norm is 
achieved as a maximum, is as follows. We consider the adjoint L* of Lo (seen as an operator 
on the Hilbert space £^([0, oo), R'")), and the composition M = L^Lo : £^([0, oo), M™) ^ 
£^([0, oo), R"*). The operator M is self-adjoint and (since Lo is a convolution operator with an 
£^ kernel) compact. Its spectrum consists of real and nonnegative eigenvalues, and its largest 
eigenvalue A is such that ji = v^A is the largest singular value of Lo, and equals the norm of Lo 
as an operator /^^([O, oo), W') -^^^([0, oo), W). Take any eigenvector u corresponding to A, 
so Mu = \u. It follows that ||Lom||^ = (LqM, Lqm) = (m, Mu) = {u, Am) = /x^ \\u\f, so m is a 
maximizing vector for Lg. Moreover, for a compact positive operator M on a Hilbert space, the 
Krein-Rutman Theorem says that, provided that there is a nonzero eigenvalue (which there is 
in this case, since M is self-adjoint and we may assume without loss of generality that M ^ 0), 
then the maximal eigenvalue A admits a nonnegative eigenvector u. Thus ||Lom|| is maximized 
at this M e i:2([0,oo),M^o)- 

6 Cascades 

Signaling systems are often built by cascading subsystems, so it is interesting to verify that a 
cascade of any number of systems which satisfy our properties again has the same form. It is 
enough, by induction, to show this for two cascaded systems 



each of which satisfies our assumptions, under the series connection obtained by setting u = 
V. The composite system can be represented in terms of the following A{^, () and B{^, Q 




x + y 



X 



z 



Ai{x)x + Bi{x)u V = Ci{x)x 
A2{z)z + B2{z)u y = C2{z)z 
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matrices: 

A-( MO 

and output y. 

It is easy to verify all the necessary properties. For example, the only nontrivial part of 
the maximization property amounts to checking that i?2(C)C'i(0 — -^2(0)6*1(0), which fol- 
lows from i?2(C)C'i(0 — -S2(0)Ci(^) (using the maximization property for B2 and the pos- 
itivity of Ci) and i?2(0)Ci(^) < i?2(0)Ci(0) (using maximization for Ci and positivity of 
i?2(0)). Similarly, the only nontrivial part of the positivity property involves checking that 
(-B2(C)C'i(0)ii ^ provided that Q = 0, for all j. But, for such a vector C, we know that 
B2iOik > for all k, so indeed B2iC)ikCiiOkj > 0. 

7 Remarks and Conclusions 

We provided a way to compute, for systems of a special form, the induced norm of the 
system. The special form includes a variety of cellular signaling cascade systems. An even 
wider class of systems can be included as well, provided that one extend our treatment to 
systems that are monotone with respect to orders other than that given by the first quadrant. 
Such orders have proven useful in analyzing, for example, MAPK cascades, see for example lEl 
m. The details of this extension will be provided elsewhere. 
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